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ABSTRACT 

We present new Gemini Planet Imager observations of the young exoplanet 51 Eridani b which 
provide further evidence that the companion is physically associated with 51 Eridani. Combining 
this new astrometric measurement with those reported in the literature, we significantly reduce the 
posterior probability that 51 Eridani b is an unbound foreground or background T-dwarf in a chance 
alignment with 51 Eridani to 2 x 10“^, an order of magnitude lower than previously reported. If 
51 Eridani b is indeed a bound object, then we have detected orbital motion of the planet between 
the discovery epoch and the latest epoch. By implementing a computationally efficient Monte Carlo 
technique, preliminary constraints are placed on the orbital parameters of the system. The current set 
of astrometric measurements suggest an orbital semimajor axis of 141^3 AU, corresponding to a period 
of 411^^2 years (assuming a mass of 1.75 Mq for the central star), and an inclination of 138l]^3 
The remaining orbital elements are only marginally constrained by the current measurements. These 
preliminary values suggest an orbit which does not share the same inclination as the orbit of the 
distant M-dwarf binary, GJ 3305, which is a wide physically bound companion to 51 Eridani. 

Subject headings: planets and satellites: detection — stars: individual (51 Eri) — planetary systems 
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1. INTRODUCTION 

Monitoring the orbi tal motion of exoplanets, through 
direc t imaging fe.g.. jChauyin et alJ 120121: iKalas et alJ 
l2013t iBonnefov et al.l 12014 iNielsen et alJ 120141) , or 
through indirect techniques su ch as radial velocity 
and transit measur e ments Ce.g.. gumming et al.l 120081: 
Howard et alJ 120121: iMarcv et alJ 12014 iMoutou et al l 


201,^ . can provide a wealth of information about their 


properties, the processes through which they form, and 
how they interact dynamically with other bodies in the 
system. Accurately determining the orbital parame¬ 
ters of exo planets can constrain the ir masses and den¬ 
sities le.g.. iCharbonneau et al.ll2000H and lead either to 
the discovery of add itional planets in the system (e.g., 
iNesvornv et al.ll2012D . or to the exclusion of additional 
planets within a range of per iods by invoking dy nami- 
cal stability arguments (e.g., iCorreia et al.l I20d5li . Or¬ 
bital parameters also provide insight as to how plane¬ 
tary companio ns dynamically interact with c ircumstellar 
material le.g.. iMillar-Blanchaer et al.ll201!^ . While the 
orbital periods are typically decades or longer for directly 
imaged planets, accurate astrometric monitoring of these 
systems can lead to preliminary constraints on their or¬ 
bital parameters before a significan t portion of the orbit 
is observed fe.g.. lPuevo et al.ll2015ll . 

51 Eridani (51 E ri) is a nearby (29.43 ± 0.29 pc; 
Ivan Leeuwen| __[2007 1 mem ber of the young (24 ± 
3 Myr; Bell et al. I2015h (3 Pictoris moving group 
( Zuckerman et al.l l2001fl . Recently, iMacintosh et al.! 
( 201511 reported the discovery of a low-mass (2-10 Mjup) 
planet at a projected separation of 13.2 ± 0.2 AU based 
on observations with the Gemini Planet Imager (GPI; 
IMacintosh et al.l 1201411 . Based on J- and i7-band spec¬ 
troscopy, and L' photometry, 51 Eri b was shown to 
have a spectrum with strong methane a nd water absorp¬ 
tion, with a temperature of 600-750 K (IMacintosh et al.l 
1201511 . Due to the short baseline between discovery 
and follow-up, it was only possible to rule out a sta¬ 
tionary background object from a 2003 non-detection; 
closer brown dwarf interlopers with non-zero proper mo¬ 
tions could not be excluded. A sta tistical argument 
based on the space densit y of T-dwarfs (iRevle et al.ll201(il : 
iBurningham et ahl 1201311 . and the allowed range of dis¬ 
tances of a foreground or background object based on the 
apparent magnitud e of 51 Eri b combined w ith the lumi¬ 
nosity of T-dwarfs (|KirkDatrick et al.1l2012[l . was used to 
limit the possibility of an unassociated field brown dwarf 
to a probability of 2.4 x 10“®. 

In addition to the newly resolved planetary companion, 
51 Eri has a n infrared excess indicative of a circumstellar 
debri s disk (|Patel et al.l 12014 iRiviere-Marichalar et al.l 
1201411 . The debris disk has yet to be spatially resolved, 
so its geometry is unconstrained. At a projected sep¬ 
aration of 1960 AU lies GJ 3305—an M -dwarf binary 
with a semimajor axis of 9.80 ± 0.15 AU (jMontet et al.l 
l2015fl —which is co-m oving with 51 Eri, for ming a bound 
hierarchical system (jFeigelson et al.l 120061) . While the 
inclination of th e GJ 3305 binary is well constrained 
{i = 92? 1 ± 0?2; iMontet et al.ll201^ . the period of the 
wide 51 Eri-GJ 3305 binary (^ 10"* years: Ipbigelson et all 
120061) precludes any such estimation of its inclination. 


2. OBSERVATIONS AND DATA REDUCTION 



Figure 1. Astrometry of 51 Eri b from 2014 December to 2015 
September (filled red, blue, and green circles). The non-moving 
background object hypothesis (light gray track), computed from 
the iifipparcos-measured parallax and proper motion of 51 Eri, is 
robustly rejected using these data alone. The measured displace¬ 
ment is well within the range of orbital motion expected for bound 
planetary-mass companions (blue envelopes). Likely orbital tracks 
were generated using the Monte Carlo method described in Sec¬ 
tion 0[to produce 10^ orbits fit to the first epoch, with the plotted 
ranges encompassing 68%, 95%, and 99.7% of the orbits. 


51 Eri was initially observed with GPI at Gemini South 
as a part of the GPI Exoplanet Survey (GPIES) on 2014 
December 18 UT (GS-2014B-Q-500). A faint companion 
candidate was identified, and subsequent observations 
demonstrated that its spectral energy distribution was 
consistent with that of a low-temperature, low-surface 
gravity giant planet (|Macintosh et al.l 1201511 . In total, 
51 Eri b has been successfully observed four times since 
discovery in 2014, three times with GPI (GS-2014B-Q- 
501, GS-2015A-Q-501), and once with NIRC2 at the 
W. M. Keck 2 telescope using the fa cility adaptive op¬ 
tics system (|Wizinowich et aP 120001) . A summary of 
these observations is given in Table [TJ For each GPI 
epoch, the observing strategy was the same. The target 
was acquired befor e transit to maximize field rotation 
(|Marois et al.ll2006l) , and the observations were taken us¬ 
ing the spectral coronagraphic mode, with either the J- 
or iJ-band Hlters. In addition to these successful obser¬ 
vations, 51 Eri was also observed with GPI in J-band on 
2015 January 29, where 51 Eri b was not recovered due 
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to poor image quality, and in i/-band on 2015 August 
30 and 31 where, although 51 Eri b was recovered, the 
signal-to-noise ratio was significantly worse. 

The GPl observations obtained on 2015 September 1 
were r educed using the GPI Data Reduction Pipeline 
(DRP; IPerrin et al.ll20llf*1 . The dark current was sub¬ 
tracted and bad pixels were identified and fixed. The 
shift in the position of the micro spectra on the detec¬ 
tor due to mechanical flexure was measured by compar¬ 
ing reference argon arcs taken monthly, to arcs obtained 
after target acquisition (|Wolff et al.11201'3) . The micro 
spectra were extracted, converting the 2D image into a 
3D {x, y, A) datacube. These were then divided by a flat 
held to correct for lenslet throughput, and were interpo¬ 
lated along the wavelength axis to a common wavelength 
vector across the bandpass. Finally, the optical distor¬ 
tion was correct ed for using measuremen ts obtained with 
a pinhole mask (iKononackv et al]l2014ll . 

To minimi ze potential biases betw een the astrometry 
presented in [Macintosh et al.l (l2015f) and the new mea- 
surements presented he re, we used Pipelines 1 and 3 from 
[Macintosh et al.l (|2015ll to both perform the point-spread 
function (PSF) subtraction, and to extract the astrome¬ 
try of 51 Eri b . Pipeline 2 was s witched to a Python im¬ 
plementation (jWang et al.l[2015f) of the Karhunen- Loeve 
Image Projection algorithm (|Soummer et al.l[2012[l . and 
uses a forward-modeled PSF to perform a Markov Chain 
Monte Carlo (MCMC) analysis to determine the poste¬ 
rior distributions of the separation and position angle of 
51 Eri b. 

The plate scale and position angle of GPI have been 
monitored by continually observing a set of astromet¬ 
ric calibrators with well-determined orbital solutions or 
contemporaneous NIRC2 me asurements, which has an 
accurate astrometric solution ([Yelda et al.|[201CI[) In ad - 
dition to observations listed in [Konopackv et ^l1 (|2014f) . 
we have observed the 9^ Ori B quadruple system an ad¬ 
ditional four times, the HD 157516 binary twice, and 
the HIP 80628 binary once. These observations were re¬ 
duced as above, and PSF subtraction was not required. 
T he pixel positions of each component were measured as 
in [Konopackv et ^ (|2014ll . 

Combining these measurements results in a plate scale 
of 14.166 ±0.007 mas lenslet”^, and a position angle off¬ 
set of —1?10 ± 0?13. The position angle offset is defined 
as the angle between the lenslet y-axis, and the elevation 
axis of the telescope, measured east from north. Since 
version 1.2, the DRP takes into account an offset of —1° 
while processing the data, and as such the difference be¬ 
tween the true position angle and the measured posi¬ 
tion angle in a DRP-reduced image is 0true — ^^measured = 
—0?10±0?13. There is no evidence of variations of either 
the plate scale or position angle offset between observing 
runs within the measurement uncertainties, so a single 
value is adopted for all epochs. The revised astrometric 
calibration was used to recalculate the astrometry from 
[Macintosh et al.l ([2015f) . which are consistent with the 
previous values, and are shown in Table [TJ 

3. COMMON PROPER MOTION CONFIRMATION 

Typically, confirmation of common proper motion is 
achieved by comparing the motion of a candidate with re- 

^ http://docs.planetimager.org/pipeline/ 
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Figure 2. Normalized posterior distributions of the proper mo¬ 
tion and distance of 51 Eri b assu ming it is an un b ound brown 
dwarf, derived from the epochs in [Macintosh et all 1)20151) using 
the updated astrometric calibration and the 2003 non-detection 
(blue histogram). The addition of the 2015 September measure¬ 
ment significantly improves the constraint on each parameter (red 
histogram). The proper motion and distance to 51 Eri are de¬ 
noted by the vertical lines. The distance posterior distribution is 
strongly peaked at the distance of 51 Eri. Using this constraint, 
the posterior probability that 51 Eri b is an unbound field dwarf is 
calculated as 2 x 10“^. 

spect to a background track for a stationary background 
object with negligible parallax, so that the only relevant 
movement is th e parallax and prop er motion of the pri¬ 
mary star (e.g.. [Nielsen et aL|[2013li . Such an analysis for 
51 Eri b is shown in Figure [TJ with the astrometry diverg¬ 
ing from the stationary background object track (2.8-(t 
in p, 8.8-(t in 9). The spectrum of 51 Eri b excluded the 
distant background star hypothesis, leaving two plausible 
hts to the spectrum: an unbou nd field brown dwarf, o r 
a bound planet, as described bv [Macintosh et aki ([201511 . 
who found a probability of 2.4 x 10“® that the object was 
an unbound field brown dwarf. By including the latest 
epoch, this probability can be further reduced by solv¬ 
ing for the allowable parallax and proper motion of an 
unbound object, and reducing the volume in which an 
unbound T-dwarf could exist. 

Using the measured par allax and proper motion of 
51 Eri (jvan LeeuwenI [2007( 1 the relative astrometry of 

51 Eri b was converted into absolute astr ometr y. _ A 

Metropolis-Hastings MCMC technique (e.g., lhbTal200ll 
was then used to fit the proper motion and parallax of 
51 Eri b to both the absolute a strometry, and the non - 
detection from 2003 discussed in [Macintosh et aki (j2015ll . 
A uniform prior in the two proper motion directions was 
used, and a p{d) oc (P prior for the distance, with a 
maximum distance of 200 pc. Good convergence was 
achieved in the MCMC chains, with a Gelman-Rubin 
statistic < 1.00007. The posteriors from the MCMC fit 
are shown in Figure |2| The median distance from the 
posteriors is 31 pc, with a 68% confidence interval be¬ 
tween 25 and 39 pc, consistent with the distance to 51 Eri 
(29.4±0.3 pc). The offset in proper motion with respect 
to the star is consistent with orbital motion occurring 
over the span of the observations. 

The calc ulation of the unbound brown dwarf proba¬ 
bility from [Macintosh et al.[ (|2015[ 1 is updated using this 
new distance constraint. The previous calculation found 
the product of the number density of T-dwarfs and the 
volume of a cone with the angular width of the GPI de¬ 
tector out to the largest distance a T-dwarf could be seen 
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Figure 3. (left): One hundred randomly selected orbits for 51 Eri b from the Monte Carlo technique described in Section |4] consistent 
with the astrometry. The color corresponds to the epoch of a given location along each orbital track, epochs later than 2055 are plotted 
in black. The position of 51 Eri is indicated by the black star. The region plotted in the right hand side is indicated by the black square 
(right): As the left panel, but focusing on the available astrometry. The color on the orbital tracks has been removed for clarity. The color 
of the symbols are as in Figure 1, and the multiple GPI ff-band epochs are labelled. 


in the reduced GPI image of a given target, repeated for 
each of the 44 stars observed as a part of GPIES at the 
time of the discovery. There was a further correction 
factor based on the proper motion constraints from the 
2003 non-detection, as only 34% of the expected distri¬ 
bution of field brown dwarfs had the motion required to 
place the object behind the star in 2003. This correction 
is not used here as the revised proper motion estimate 
for the unbound brown dwarf scenario is S-cr discrepant 
from the proper motion of the star due to orbital mo¬ 
tion (Figured]). Instead, the volume calculation of the 
cones was convolved with the distance posterior distri¬ 
bution, lowering the posterior probability that 51 Eri b 
is an unbound brown dwarf to 2 x 10“^. 

4. ORBITAL ELEMENTS OF 51 ERI b 

MGMG orbit fitting is slow to converge for sparsely 
sampled astrometry or short orbital arcs, and so we im¬ 
plemented a more computationally efficient Monte Carlo 
method to generate plausible distributions of orbital pa¬ 
rameters based on astrometry covering only a small frac¬ 
tion of an orbital period. For four of the orbital pa¬ 
rameters, a large number of samples were drawn from 
appropriate probability distributions: uniform for argu¬ 
ment of periastron (ui) and epoch of periastron passage 
(To), uniform in cos(z) for inclination angle (i), and ec¬ 
centricity (e) following t he lin ear fit to radial velocity 
planets of iNielsen et al.l (I2008D . Initial orbits were gen¬ 
erated using these four parameter distributions and fixed 
values of semimajor axis (a) and position angle of nodes 
(n), the values of a was scaled and the value of 12 rotated 
to reproduce the astrometry at the first epoch. Period 
(P) was not fit a nd was calculated assu ming a stellar 
mass of 1.75 Mq (iSimon fc Schaefeill2011l) . Astrometric 
errors were incorporated into the generated parameters 


by adding random offsets to the separation and position 
angle of the first epoch before each orbit is shifted and 
rotated. These uncertainties were randomly drawn from 
Gaussian distributions equivalent to the first epoch as¬ 
trometric uncertainties. 

With multiple epochs, we proceeded to iteratively re¬ 
ject sets of orbital parameters that do not match later 
measurements, with acceptance probability given by a 
two-dimensional Gaussian with the astrometric errors as 
the standard deviations. Generated orbits that were clos¬ 
est to the observed separation and position angle at the 
corresponding observational epoch were more likely to be 
accepted. We then obtained distributions of fitting or¬ 
bital parameters given the input astrometry. Varying the 
epoch chosen to initialize the procedure had little effect 
on the distribution of accepted orbital parameters. 

We validated this method by generating 10^ orbits with 
one fixed orbital parameter and randomly sampling the 
other parameters, and then creating five epochs from one 
year of simulated astrometry for each orbit, each with an 
observational uncertainty. The spacing of the epochs and 
astrometric errors were chosen to be the same as the mea¬ 
surements of 51 Eri b. We then applied our method to the 
artificial astrometry and examined the returned distribu¬ 
tions of orbital parameters. Of the simulated orbits with 
semimajor axis of 13 AU, the median was 13.3 AU, and 
in 74% of orbits 13 AU was within the 68% confidence 
interval of each individual orbital fit for semimajor axis. 
This suggests the generated distributions are reasonable 
representations of the posterior probability distributions, 
and comparison with MGMG (described below) further 
supports this. This procedure will be described further 
in Blunt et al. (2015, in preparation). 

This technique was applied to the astrometry of 
51 Eri b. One hundred orbits from the fit are shown in 
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Figure 4. Normalized distributions of accepted orbital parameters consistent with the astrometry from the Monte Carlo a nalysis described 
in Sec tion m The analysis was performed using the discovery epoch (black dashed histogram), using all epochs presented in IMacintosh et ahl 
(green dotted histogram), and using all epochs (red solid histogram). As more astrometry is added, the inclination angle becomes 
significantly more constrained. The normalized posterior distributions from the MCMC analysis are also shown for comparison (blue 
dotted-dashed histogram). 


FigurelSl and the distributions of accepted orbital param¬ 
eters are shown in Figured While the marginal distribu¬ 
tion for inclination using the discovery epoch follows the 
shape of the prior, the distribution after incorporating 
data from 2015 peaks at ^ 138°, showing the new ob¬ 
servations have provided significant orbital constraints. 
More moderate changes are seen in eccentricity and semi¬ 
major axis, with less eccentric orbits, and orbits closer to 
14 AU, becoming more favored with increasing number 
of measurements. The orbital elements corresponding to 
the minimum , representing the best-fitting orbit gen¬ 
erated, are given in Table [TJ along with the median and 
limits corresponding to where 68% of the Monte Carlo 
generated orbits were found. In order to compute the 
median and 68% intervals, ui was wrapped to be within 
0°-180°, n to within 30°-120°, and Tq to within 1995- 
1995-bP. 


The shapes of the distributions of accepted orbital pa¬ 
rameters were confi rmed using an affine invariant MCMC 
ensemble sampler (|Foreman-Mackev et al.l[2013fFI . with 
the same prior on each parameter. The shapes of the 
posterior distributions are similar (Figure 0]); the me¬ 
dian and l-cr confidence intervals for a and i were found 
to be Id.lj'jg'Q AU and 135lg3 consistent with the 
values in Table [U but taking two orders of magnitude 
longer to compute. In addition to these Monte Carlo 
techniques, the astrometry of 51 Eri b was applied to 
the technique for constrainin g orbital parameter s over 
short orbital arcs presented in lPearce et al.l (l2015f ). The 
angle between the projected separation and velocity vec¬ 
tors was calculated as ^ = 84 . 5 ^ 14(3 deg, and a value 
for their dimensionless parameter B of 0.23lgjg. Com- 

^ http://dan.iel.fm/emcee 
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Table 1 

Observations and Orbital Parameters of 51 Eridani b 


UT Date 

MJD 

Instrument 

Filter tint 

-^coadd 

A^exp 

Plate Scale 

Position Angle 

P 

e 




(s) 



(mas px~^)^ 

Offset (deg) 

(mas) 

(deg) 

2014 Dec 18 

57009.1292 

GPI 

H 59.6 

1 

38 

14.166 ±0.007 

-1.10 ±0.13 

449.8 ±6.7 

171.0 ±0.9 

2015 Jan 30 

57052.0572 

GPI 

J 59.6 

1 

70 

14.166 ±0.007 

-1.10 ±0.13 

454.0 ± 6.4 

170.6 ± 1.0 

2015 Jan 31 

57052.9753 

GPI 

H 59.6 

1 

64 

14.166 ±0.007 

-1.10 ±0.13 

461.8 ± 7.1 

170.5 ±0.9 

2015 Feb 01 

57054.0364 

NIRG2 

L' 0.9 

60 

62 

9.952 ± 0.002‘> 

0.252 ± O.OOO*’ 

461.5 ± 23.9 

170.4 ±3.0 

2015 Sep 01 

57266.4052 

GPI 

H 59.6 

1 

93 

14.166 ±0.007 

-1.10 ±0.13 

454.7 ±5.7 

166.5 ±0.6 




Preliminary Orbital Parameters ot 51 Eridani b 




Parameter 

Unit 

■ 

Median 

Lower 

Upper 



Semimajor 

axis ia) 

AU 

18.61 

14 


11 

21 



Eccentricity (e) 

- 

0.2804 

0.21 


0.06 

0.40 



Inclination (i) 

deg 

126.7 

138 


125 

153 



Argument of periastron (a;) 

deg 

4.233 

90 


32 

148 



Position angle of nodes (D) 

deg 

166.6 

75 


47 

100 



Epoch of periastron (To) 

- 

2016.05 

2014.79 

2003.31 

2026.11 



Period (P) 

year 

60.68 

41 


29 

76 




^ lYelda et'HI (f20Toh 

In reduced GPI datacubes, one pixel is equivalent to one lenslet 

paring these to the minimum inclin ation and eccentric- 
ity contours of iPearce et al.l (j2015fl . the eccentricity of 
51 Eri b is unconstrained, and the inclination is restricted 
to i < 771y deg (or i > lOStg deg). While consistent 
with the values in Table [H these limits are significantly 
less constraining. 

5. CONCLUSIONS 

We have presented a new astrometric measurement of 
the young exoplanet 51 Eridani b, which provides fur¬ 
ther evidence supporting the bound companion hypoth¬ 
esis. Using this new measurement, we observe significant 
deviation {2.8-a in p, 8.8 -(t in 9) from the stationary 
background object track, with a reduced Xhkg 11-19 
(|Nielsen et al.ll20l'll . We compute a probability that it is 
an interloping field brown d warf of 2 x 10~^, d e crease d by 
a factor of twelve relative to [Macintosh et al.l (|2015f) . By 
implementing a computationally efficient Monte Carlo 
method to sample probable orbits, we place the first con¬ 
straints on the orbital parameters of 51 Eri b. Based 
on the present astrometry, the median of the semima¬ 
jor axis distribution is 14 AU, corresponding to a period 
of 41 years, assuming a mass of 1.75 Mq for 51 Eri. 
While the additional astrometric epoch did not signif¬ 
icantly change the distribution of accepted semimajor 
axes and periods, the range of allowed inclinations was 
significantly constrained, with a median value of 138°. 
The eccentricity of the orbit remains unconstrained, with 
circular orbits only marginally preferred relative to the 
prior distribution. 

Based on these preliminary constraints, the orbit of 
51 Eri b does not appear to be co-planar with the or¬ 
bit of the M-dwarf binary GJ 3305 (i = 92? 1 ± 0?2; 
iMontet et al.l[20T5ll . The inclination of the orbital plane 
of the outer binary relative to 51 Eri is not known. 
The large separation between 51 Eri and GJ 3305, and 
the young age of the system, would suggest that sec¬ 
ular Lidov-Kozai oscillations would not have had suf¬ 
ficient time t o significantly alter the semimajor axis 
of the planet (IMontet et al.l [2M^ . although moderate 
changes in the inclination and eccen tricity are not ex¬ 
cluded (jEabrvckv fc Tremainel 1200711 . Given the long 


timescale of these oscillations (2(10 My r for a perturber 
on a circular orbit. IMontet et aLll2015ll . this effect could 
be completely suppressed by quicker secular precession 
due to, e.g., a relat ively low-mass, undetected planet 
(IWu fc Murravl[200^ . Alternatively, the planet may re¬ 
main o n its current orbit, despite the presence of the wide 
binary (jHolman fc Wiegertill999ll . Gontinued astromet¬ 
ric monitoring of 51 Eri b over the next few years should 
be sufficient to detect curvature in the orbit, further con¬ 
straining the semimajor axis and inclination of the orbit, 
and placing the first constraints on the eccentricity. Ab¬ 
solute astrometric measure ments of 51 Eri with GAIA 
(e.g., iPerrvman et al.l 1201411 . in conjunction with moni¬ 
toring of the relative astrometry of 51 Eri b, will enable a 
direct measurement of the mass of the planet. Gombined 
with the well-constrained age of 51 Eri b, such a deter¬ 
mination would provide insight into the evolutionary his¬ 
tory of low-mass directly imaged extrasolar planets, and 
help distinguish between a hot-start or core accretion 
formation process for this planet. 


Based on observations obtained at the Gemini Ob¬ 
servatory, which is operated by the Association of Uni¬ 
versities for Research in Astronomy, Inc., under a co¬ 
operative agreement with the National Science Eounda- 
tion (NSF) on behalf of the Gemini partnership: the 
NSF (United States), the National Research Gouncil 
(Canada), CONICYT (Chile), the Australian Research 
Council (Australia), Ministerio da Ciencia, Tecnologia 
e Inovagao (Brazil) and Ministerio de Ciencia, Tec- 
nologfa e Innovacion Productiva (Argentina). This re¬ 
search has made use of the SIMBAD database, op¬ 
erated at CDS, Strasbourg, France. Supported by 
NSF grants AST-0909188 and AST-1313718 (R.J.D.R., 

J. R.G., J.J.W., T.M.E., P.G.K.), AST-1411868 (B.M., 

K. F., J.L.P., A.R., K.W.D.), AST-141378 (P.A., 

G.D., M.P.F.), NNX11AF74G (A.Z.G., A.S.), and 
DGE-1232825 (A.Z.G.). Supported by NASA grants 
NNX15AD95G/NEXSS and NNX11AD21G (R.J.D.R., 
J.R.G., J.J.W., T.M.E., P.G.K.), and NNX14AJ80G 
(E.L.N., S.C.B., B.M., F.M., M.P.). J.R., R.D. and 
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